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Percolation thresholds of the Fortuin-Kasteleyn cluster for the Edwards-Anderson 
model on complex network: analytical results on the Nishimori Line 

Chiaki Yamaguchi 

Kosugichou 1-359, Nakahara-ku, Kawasaki-si, Kanagawa 211-0063, Japan 

We show analytically the percolation thresholds of the Fortuin-Kasteleyn cluster for the 
Edwards- Anderson model on random graphs with arbitary degree distributions. The results 
on the Nishimori line are shown. The results for the ±J model, the diluted ±J model, and 
the Gaussian model are respectively shown. The results are shown by applying an extension 
of a criterion for the random graphs with arbitary degree distributions. The results for the 
infinite-range ±J model and the Sherrington-Kirkpatrick model are also shown. 
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1. Introduction 

Study of spin models on complex network has been 
important 1 with the development of study of complex 
network. In the present article, random graphs with ar- 
bitary degree distributions are adopted as examples of 
the complex network. In the present article, behavior of 
spins on no growing network is investigated. 

We investigate the Edwards-Anderson model 2 as an 
Ising spin-glass model. The understandings for the 
Edwards- Anderson model on random graphs and on the 
Bethe lattice are still uncompleted. 1 ' 3,4 In the present ar- 
ticle, the ±J model, the diluted ± J model, and the Gaus- 
sian model in the Edwards- Anderson model are respec- 
tively investigated. The Nishimori line in the Edwards- 
Anderson model is a line in the phase diagram for the 
exchange interactions and the temperature. The inter- 
nal energy, the upper bound of the specific heat, and so 
forth are exactly calculated on the Nishimori line in the 
Edwards- Anderson model. 5-9 The location of the mul- 
ticritical point in the Edwards- Anderson model on the 
square lattice is conjectured, and it is shown that the 
conjectured value is in good agreement with the other 
numerical estimates. 10 In the present article, the results 
on the Nishimori line are shown. 

There is a case for occuring a percolation transition 
of network that the network is devided into a lot of net- 
works by deleting nodes and/or links on a network. We 
define the percolation transition as the percolation tran- 
sition of network in the present article. There is a case for 
occuring a percolation transition of cluster that the clus- 
ter consists of bonds put between spins becomes a giant 
component. The bond is fictitious. We define the percola- 
tion transition as the percolation transition of cluster in 
the present article. In the present article, the percolation 
transition of cluster on a complex network is mentioned. 

In the present article, the percolation transition of the 
Fortuin-Kasteleyn (FK) cluster in the FK random cluster 
model 11 is investigated. In the ferromagnetic spin model, 
the percolation transition point of the FK cluster agrees 
with the phase transition point. The agreement gives the 
geometrical understanding of the phase transition phe- 
nomenon. 12 Powerful Monte Carlo methods using the FK 



cluster have been proposed. 13-17 On the other hand, in 
the Edwards- Anderson model that has a conflict in the 
interactions, the percolation transition point of the FK 
cluster disagrees with the phase transition point. The dis- 
agreement makes the geometrical understanding of the 
phase transition phenomenon awkward. There are plenty 
of approches for solving the awkwardness by using exten- 
sions of the FK random cluster model. 18 On the other 
hand, it is pointed out, by Arcangelis, A. Coniglio, and 
F. Peruggi, that the correct understanding of the per- 
colation phenomenon of the FK cluster in the Edwards- 
Anderson model is important since a dynamical transi- 
tion, which is characterized by a parameter called the 
Hamming distance or damage, is occurred at a temper- 
ature very close to the percolation temperature, and the 
dynamical transition and the percolation transition are 
related to a transition for a signal propagating between 
spins. 19 In the present article, the percolation threshold 
of the FK cluster in the Edwards- Anderson model is an- 
alytically shown. 

Study of the random graphs with arbitary degree dis- 
tributions has been developed. 20 Our results are shown 
by applying an extension of a criterion 2 1-23 for the ran- 
dom graphs with arbitary degree distributions. The re- 
sults for the infinite-range ±J model and the results 
for the Sherrington-Kirkpatrick (SK) model 24 are also 
shown. 

The present article is organised as follows. First, 
a complex network model and the Edwards- Anderson 
model are described. The FK cluster is described. A cri- 
terion for percolation of cluster is described. The per- 
colation thresholds are shown for the ±J model and the 
diluted ± J model. The percolation threshold is shown for 
the Gaussian model. The present article is summarized. 

2. A complex network model and the Edwards- 
Anderson model 

Network consists of nodes and links connected between 
nodes. In the present article, the network models that 
we investigate are random graphs with arbitary degree 
distributions. The network has no correlation between 
nodes. The node degree, k, is generated according to a 
distribution of the node degree, p{k). The links are ran- 
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Fig. 1. A relation between the aspect a of the network and the 
model on the corresponding network of the aspect a. 



domly connected between nodes. 

We define a variable b(i, j) where b(i, j) gives one when 
node i and node j is connected by a link and gives zero 
when node i and node j is not connected by the link. 
The degree k(i) of node i is given by 



(1) 



The average of the node degree for links, < k >jy, is 
given by 



1 N 



(2) 



where < > n is the average over the whole network. N is 
the number of nodes. The average of the square of the 
node degree for links, < k 2 >n, is given by 



1 N 

< fc2 >^E fc2 (*)- 



We define 



2<k> 



N 



(3) 



(4) 



< k 2 > N 

where a represents an aspect of the network. 

Fig. 1 shows a relation between the aspect a of the 
network and the model on the corresponding network of 
the aspect a. A network is almost the complete graph 
when a is close to zero, and the model on the network 
is almost an infinite-range model. The model on the net- 
work is the infinite-range model when < k >n= N — 1, 

< k 2 > N = (N - l) 2 , and a = 2/(N - 1). We define 
a coordination number as z. p(k) = 5k, z where S is 
the Dirac delta function. The coordination number z is 
two when the aspect a of the network is one. The net- 
work consists of many cycle graphs when z is two. The 
model on the network consists of many chain models. 
In the Erdos-Rcnyi (ER) random graph model and in 
the Gilbert model, the distribution of node degree is the 
Poisson distribution. 1 The ER random graph model is 
a network model that the network consists of the fixed 
number of nodes and the fixed number of links, and the 
links are randomly connected between the nodes. The 
Gilber model is a network model that link between nodes 
is connected with a given probability. In the ER random 
graph model and in the Gilbert model, < k >n~ 1 and 

< k 2 >at=< k >jy (< k >n +1) = 2 when a is one. 
The Hamiltonian for the Edwards- Anderson model, H, 



is given by 



N 



n 



1 {j\b(i,j)=l} 



(5) 



where Sj = ±1. Jij is the strength of the exchange inter- 
action between the spin on node i and the spin on node 
j. The value of Jij is given with the distribution P(Jij). 
The ±J model, the diluted ±J model, and the Gaus- 
sian model are respectively given with the difference of 

P(Jij)- 

For the ±J model, the distribution P( ±J \Jij) of Jij 
is given by 



p( ±J Xj ij )=p6j ij> j + (l-p)6 Jij ,-j 



(6) 



where J > 0. p is the probability that the interaction 
is the ferromagnetic interaction (J^ = J). 1—p is the 
probability that the interaction is the antiferromagnetic 
interaction (Jij = — J). 

For the diluted ±J model, the distribution 
P( D±J )(J ij ) of Jij is given by 

P (U±J) (Jij) +qS Jij ,-j +r5j ij , (7) 

where J > and p + q + r = 1. p is the probability that 
the interaction is the ferromagnetic interaction ( Jij = J). 
q is the probability that the interaction is the antiferro- 
magnetic interaction ( = — J). r is the probability that 
the interaction is diluted (Jij = 0). This model is the ± J 
model when r = 0. 

For the Gaussian model, the distribution 

p(Gaussian)(j.^ Q f j.. ig giyen by 

p(Gaussian) ^ j. ^ _ \ Q - (Jg - J ) 2 /2.7 2 ^ ^g-j 

V2tt J 2 

The average of Jij is {Jij}n = Jo where []r is the 
random configuration average. The variance of Jjj is 

[ J ij]R - [Jii?R = J 2 - 

In the present article, for calculating thermodynamic 
quantities, a gauge transformation: 



Jij ^ J ij @i @ j 



Si — > Si* 



(9) 



is used where Oi = ±1. Using the gauge transformation, 
7i — > H and P(Jij) — > P(Jij<Ji<jj). This transformation 
has no effect on thermodynamic quantities. 25 

3. The Fortuin-Kasteleyn cluster 

The bond for the FK cluster is put between spins 
with the corresponding probability Pbond(5'i, Sj , Jij) of 
the spin states and the strength of the exchange interac- 
tion. The probability Pbond(Si> Sj, Jij) is given by 

Pbond(Si, Sj, Jij) = i - e-ew'-ftw (10) 

19 where /3 is the inverse temperature and (3 — 1/ksT. 
ks is the Boltzmann constant and T is the temperature. 
By connecting the bonds, the FK clusters are generated. 

The thermodynamic quantity of the bond put be- 
tween the spin on node i and the spin on node j, 
[< b hond (i, j) >t]r, is 

[< b hond (i, j) > T }r = [< Pbond(Si, Sj, Jij) > t }r (11) 

where < >t is the thermal average. The thermody- 
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namic quantity of the node degree for bonds at node 

i, [< fcbond(«) >t]r, is 

[< k hond (i) > T ] H = 

[< E M^Sj.JtflHfl. (12) 

{j|6(i,j)=l} 

The thermodynamic quantity of the square of the node 
degree for bonds at node i, [< k^ ond (i) >t]r, is 

[< fcbondM >T}R = 

[< y ' y ' P~bond(Si, Sj, Jij) X 

{j\b(i,j)=l} {l\b(i,l) = l} 

Pbond{S t , Si, J,()(l - + 

E f\>ond(Si, Sj, Jjj) >t]k ■ (13) 
{J'|6(»,J') = 1} 

The thermodynamic quantity of the node degree for 
bonds, [< fcbond >t]r, is 

1 W 

[< fcbond >T}R = fcbond(i) >t]r ■ (14) 

The thermodynamic quantity of the square of the node 
degree for bonds, [< fcg ond >t]r, is 

1 N 

[< fcbond >T]ii = fc bondW >T]i? • (15) 

i 

4. A criterion for percolation of cluster 

The percolation of the random graphs with arbitary 
degree distributions is occured when 



is 



< k ■ >n> 2 < k > 



N 



(16) 



21-23 j nC q_ (15) j s gi vcn by the inequality when the 
network is percolated. Ineq. (16) is given by the equal- 
ity when the network is at the percolation transition 
point. The criterion (ineq. (16)) is true for sufficiently 
large number of nodes. Ineq. (16) is derived, with var- 
ious methods, by the group of Molly and Reed, 21 the 
group of Cohen et. al., 22 and the group of Newman et. 
al., 23 respectively. 

From ineq. (16), the network is percolated when a < 1. 
From ineq. (16), the network is at the percolation tran- 
sition point when a = 1. The network and the cluster on 
the network are unpercolated when a > 1. Therefore, the 
percolation of cluster is investigated for < a < 1. 

In the case that links and/or nodes are randomly di- 
luted on the random graphs with arbitary degree dis- 
tributions, the criterion (ineq. (16)) is applicable to the 
diluted network. 22 In the case that links are diluted on 
the random graphs with arbitary degree distributions, 
the percolation problem for the diluted network can be 
regarded as the bond percolation problem. We define the 
bond states as the graph G. In the graph G, we define 
the node degree for bonds at node i as fcbond (G, i). Since 
the bonds are randomly put on the links that are ran- 
domly connected between the nodes, the criterion of the 
percolation of cluster for the bond percolation problem 



i i 

In what follows, a criterion of the percolation of cluster 
for spin models is given based on the above discussion. 

We consider a condition that the bond does not depend 
on the size of the degree k(i). We define a variable for 
the inverse temperature as p((3). We set 



< p{fi) < 1 . 



(18) 



We assume the case that [< &bond(*,j) >t]r, [< 
fcbond(«) >t]r, and [< fcb ond («) >t]r are respectively 
written in 



[< b hond (i,j) > t ]r = p(/3) , 
[< fcbond(i) >t]r = p(fi) k(i) , 



(19) 
(20) 



[< fcbond M >T}R 



p 2 03)fc(*)[fc(z)-l] + 

p{(3)k(i). 



(21) 

In the case, it is implied that the bias for the size of k(i) 
does not appear in the statistical results of bonds put 
between spins. Therefore, we describe the case that [< 
&bond(«, j) >t]r, [< fcbond (i) >t]r, and [< fc 2 ond (i) >t]r 
are respectively written in eqs. (19, 20, 21) as the case 
that the condition that the bond does not depend on the 
size of k(i) is satisfied. 

In the case that the condition is satisfied, as an exten- 
sion of ineq. (17), we obtain 



^E fc bo„ d ({^},{J,/},G,*)> 



N 



E fc bond({^},{J, ; },G,*) (22) 



since the bonds that do not depend on the size of k(i) 
are put on the links that are randomly connected be- 
tween nodes. fcbond({<Sj}, { Jji}, G, i) is the node degree 
for bonds at node i in the graph G that is compatible of 
{Sj} and {Jji}- It is anticipated that ineq. (22) is true 
for sufficiently large number of nodes in the case that 
the condition is satisfied. Using ineq. (22), we obtain the 
criterion of the percolation of cluster for spin models as 



[< fcbond >T]R > 2[< fcbond >t]r ■ 



(23) 



Ineq. (23) is given by the inequality when the cluster is 
percolated. Ineq. (23) is given by the equality when the 
cluster is at the percolation transition point. Ineq. (23) 
gives the percolation threshold of the cluster in the case 
that the condition is satisfied. 

5. The ±J model and the dilued ±J model 

The distribution p( ±J ) (J y ) of J^ for the ±J model 
is, using eq. (6), 

(±7) g/^P^ij 



2cosh(/3 P J) ' 



(24) 
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where (3p is given by 



a 1 i P 
Pp = 777 l n ■ 



2J 1 -p 



(25) 



The thermodynamic quantity of the bond put be- 
tween the spin on node i and the spin on node j, 
[< &bond >t] ( r J) , is, using eqs. (9, 10, 11, 24) when 

[< b hond (i,j) >t] ( r J) 

{,J lm } <lm> 

E{5 i} PbondjSj, Sj, Jjj)e^ J "" s ' s - 

{Si} 



2 N [2 cosh((3 P J)] N B 

]T ]T p bo„d(5 2 , Sj, Ja)^" 



= tanh(/3 P J) 



(26) 



where < xy > denotes the nearest neighbor pairs con- 
nected by links, Nb is the number of all links, and 
Nb = N < k >n /2. The thermodynamic quantity of 

the node degree for bonds at node i, [< fchond W > t]p ±J \ 
is, using eqs. (9, 10, 12, 24) when (i = ftp, 



[< k hond (i) >t] ( r J) = tanh(/3 P J) k(i) 



(27) 



The thermodynamic quantity of the square of the node 
degree for bonds at node i, [< k^ ond (i) > t]p ±J \ is, using 
eqs. (9, 10, 13, 24) when (3 = (3 P , 

[< fcbondW >t] { r J) = tanh 2 (/3 P J)fc( J )[fc( J )-l] + 

tanh^p J) k(i) . (28) 



We set 



p (±J) (/3 P ) =tanh(^pJ). (29) 



Eqs. (26, 27, 28, 29) are formulated as eqs. (19, 20, 21), 
respectively. Therefore, the bond does not depend on the 
size of k(i). Using eqs. (14, 15, 23, 27, 28), we obtain the 
criterion of the percolation of the FK cluster for the ± J 
model as 

2 < k > N 



1 - exp(-2/3 P J) > 



< k 2 > 



N 



(30) 



Ineq. (30) is given by the inequality when the cluster is 
percolated. Ineq. (30) is given by the equality when the 
cluster is at the percolation transition point. 

From eqs. (18, 29), it is realized that there is the per- 
colation transition point for < (3 P < oo. From eq. (30), 
it is realized that there is the percolation transition point 
for < a < 1. The probability p( ±J ) that the interaction 
is the ferromagnetic interaction is, using eqs. (25, 30), 



P 



(±J) 



1 



(31) 



at the percolation transition point. The percolation tran- 



sition temperature T p ±J ^ is, using eqs. (25, 31), 



T (±J) _ 
1 P — 



(32) 



In Ref., 5 the internal energy on the Nishimori line for the 
±J model, E( ±J \ is derived. The internal energy 
is equivalently 

E[±J) = jVJtanh(/3 P J) 



< k > 



N 



The internal energy is, using eqs. (32, 33), 

NaJ 



E (±J) 



< k > 



(33) 



(34) 



2(2 - a) 

We define the specific heat on the Nishimori line for the 
±J model as C {±J \ In Refs., 5 - 6 the upper bound of the 
specific heat C^ ±J ^ is derived. The upper bound of the 
specific heat C^ ±J ^ is equivalently 



C(±J) < fcpjV(/3pJ) 2 sech 2 (/3pJ) < k 



N 



(35) 



The upper bound of the specific heat C^ ±J ^ is, using 
eqs. (32, 35), 



ci±J) < k B N(l - - a) 

~ 2(2 -a) 2 V 1 



< k > 



N 



(36) 



It is realized that the specific heat has no singularity. 

In Fig. 2, the percolation threshold of the FK cluster 
for the ± J model is shown. In (a) of Fig. 2, the relation 
between the aspect a of the network and the probability 
p(± J ) that the interaction is the ferromagnetic interac- 
tion is shown. Eq. (31) is used for showing (a) of Fig. 2. 
In (b) of Fig. 2, the relation between the aspect a of 
the network and the percolation transition temperature 
T p ±JS> is shown. Eq. (32) is used for showing (b) of Fig. 2. 
In (b) of Fig. 2, J/k B = 1 is set. 

For the ferromagnetic Ising model on the same net- 
work, the phase transition temperature j^ Ferro ) [ s 



-n(Fcrro) 



fcpln^ 



(37) 



26,27 rp{±J) 



(eq. (32)) coincides with j^ Ferro ) . For the 
ferromagnetic Ising model on the same network, the spe- 
cific heat diverges or jumps at j , ( l Forro \ 2 6,27 Q n tne |; ner 
hand, from eq. (36), the specfic heat neither diverge nor 
jump at T p ±J ^ on the Nishimori line. 

The complete graph is considerd as a case that a ~ 0. 
We set < k > N = N-l, < k 2 > N = (iV-1) 2 , a = 2/(JV- 
1), and J — > J/*jN. From the settings, the model on the 
network becomes the infinit-range ±J model. In what 
follows, our results for the infinit-range ±J model are 
shown. The probability p( IR± J ~> that the interaction is the 
ferromagnetic interaction is, using eq. (31) for sufficiently 
large number of nodes, 

N-l 1 



P 



(IR±J) 



2(N - 2) 



(38) 



at the percolation transition point. The percolation tran- 
sition temperature T p m±J ^ is, using eq. (32) for suffi- 
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Fig. 2. The percolation threshold of the FK cluster for the ±J 
model, (a) The relation between the aspect a of the network and 
the probability p^ JS> that the interaction is the ferromagnetic 
interaction is shown, (b) The relation between the aspect a of 
the network and the percolation transition temperature T^ J ^ 
is shown. J/fca = 1 is set. 



ciently large number of nodes, 

rp(m±J) _ 



J_ 

k B 



N . 



(39) 



The internal energy ij( IR±J ) is, using eq. (34) for suffi- 
ciently large number of nodes, 



E (1K±J) = _ 



N(N - 



(40) 



2{N-2) ' 2 

For sufficiently large number of nodes, the right-hand 
side of eq. (36) reduces to 



k B N(N -1) 2 {N -3) 
8(iV-2) 2 



In 1 



N-3 



fee 
~~2 



The upper bound of the specific heat C^ IR±J ^ is, using 
eqs. (36, 41), 



~ 2 



(42) 



In Rcf., 18 the percolation transition temperature of the 
FK cluster for the infinit-range ± J model is derived by 
using the analytical solution of the SK model. The per- 
colation transition temperature of the FK cluster for the 
infinit-range ±J model obtained in the present article 
agrees with the result for single replica case in Ref. 18 

The case for a = 1 is described. Using eq. (31), 



p = 1. Using eq. (32), Tp = 0. From eq. (16), it is 
realized that the network is at the percolation transition 
point. From p = 1, the exchange interaction is only the 
ferromagnetic interaction. From Tp = 0, all the spins 
are parallel. Therefore, the cluster and the network are 
at the percolation transition point since the bonds are 
put on all the links of the network. 

The distribution p( D±J ) ( J y ) of J l3 for the diluted ±J 
model is, using eq. (7), 



e ( p ) J 2 +0 P J + i + e P¥>J 2 -f3 P J 



(2) 



(2) 

where f3 p ' and [3p are respectively 



(43) 



(44) 



. 9 This model becomes the ±J model when r = 0. In 
what follows, the result for r ^ is only described since 
the result for the ± J model is already desicribed above. 

The thermodynamic quantity of the bond put be- 
tween the spin on node i and the spin on node j, 
[< b hond (i,j) >t] { r ±J \ is, using eqs. (9, 10, 11, 43) when 

[< b hond (i,j) >t}r ±J) 

= E II P {B±J) (Jlm) X 
{,/,„} <lm> 

E{S i} PbondjSi, Sj, Jjj) ^<^> J '- S ' s " 



e Pp E<; m > JlmSlSr, 



1 



2A r ( e /3p ) ^ 2 +/3p./ + 1 + e ^p } J 2 -0pJjN B 

y ] y ] Pbond(Si, Sj, Jij) X 
{Jl m }{S,} 

e Pp } E< !m > JL+Pp E< !m > Ji m s l s m 



(1-r) tanh(/3 P J). 



(45) 



The thermodynamic quantity of the node degree for 
bonds at node i, [< fcbond(i) >t]p^ ±J \ is, using eqs. (9, 
10, 12, 43) when /3 = /3 P , 

[< fcbond(«) >T] { p ±J) = (1 - r) tanh(/3 P J) k(i) . (46) 
(41) The thermodynamic quantity of the square of the node 



degree for bonds at node i, [< k^ ond (i) 
using eqs. (9, 10, 13, 43) when = 0p, 



>T 



(D±J) 
R 



is, 



[< fcbondW >T] { p ±J) = 

(1 - rf tanh 2 (/3 P J) k(i)[k(i) - 1] + 
(1-r) tanh(/3 P J) k{i) . 



We set 



SD±J) 



(1 - r) tanh(/3 P J) . 



(47) 



(48) 



Eqs. (45, 46, 47, 48) arc formulated as eqs. (19, 20, 21), 
respectively. Therefore, the bond does not depend on the 



J. Phys. Soc. Jpn. 



Full 



size of k(i). Using eqs. (14, 15, 23, 46, 47), we obtain 
the criterion of the percolation of the FK cluster for the 
diluted ± J model as 

-20 P Jj 



Paper 

a). Using ineqs. (51, 57) 
1 -r 



2(1 - r)(l - e 



> 2<k>N . (49) 



< k 2 > 



N 



1 + e~ 2 Pr J + (1 - r )(l - e~ 2 ^ J ) 

Incq. (49) is given by the inequality when the cluster is 
percolated. Ineq. (49) is given by the equality when the 
cluster is at the percolation transition point. 

From eqs. (18, 48), it is realized that there is the per- 
colation transition point for r ^ 1 and < j3p < oo. 
Using incq. (49), 

(2-a)(l-r)-a 



> e 



-2f) P J 



> 0. 



(50) 



(2 - a)(l -r) + a 

Using the left-hand side of ineq. (50) and the right-hand 
side of ineq. (50), 

a 



1 - r > 



(51) 



When ineq. (51) is satisfied, there is the percolation tran- 
sition point. 

The probability p( D±J ) that the interaction is the fer- 
romagnetic interaction is, using eqs. (44, 49), 

(D±J) _ (2- a )(l-r) + a 
P ~ 2(2-o)(l-r) 1 ' 

at the percolation transition point. The percolation tran- 
sition temperature j p D±J ) is, using eqs. (44, 52), 



T (T>±J) _ _J_ f 

P k j (2-a)(l-r)+a 

111 (2-a)(\-r)-a 



(53) 



In Ref., 9 the internal energy on the Nishimori line for 
the diluted ±J model, E^- D±J \ is derived. The internal 
energy ij( D±J ) is equivalently 



E (D±J) = _ 



NJ(l-r) tanh(/3 P J) 



< k > 



N 



The internal energy i?( D±J ) is, using eqs. (53, 54), 

NaJ 



£(D±J) = _. 



< k > 



N 



(54) 



(55) 



2(2 -a) 

We define the specific heat on the Nishimori line for the 
diluted ±J model as C^ ±J \ In Ref., 9 the upper bound 
of the specific heat C' D±J ) is derived. The upper bound 
of the specific heat C^ D±J ^ is equivalently 

k B N{(3 P J) 2 



< k > N X 



[(1 - r) 2 sech 2 (/3p J) + r(l - r)} . (56) 

The upper bound of the specific heat C^ D±J ^ is, using 
eqs. (53, 56), 



C (B±J) 



< 



k B N 



(2 -a) 



In 



(2-a)(l-r)+a 



<k> N (57) 



(2-a)(l-r)-aj 

The upper bound of the specific heat diverges when r = 
2(1 - a)/ (2 - a). Using eq. (52), p=l when r = 2(1 - 
a)/(2-a). Using eq. (53), T P = when r = 2(l-a)/(2- 



2o(l-a) 
(2 - a) 2 ~ (2 -a) 2 



> 



Author Name 



(58) 



The upper bound of the specific heat C^ D±J ' is nonneg- 
ative. 

6. The Gaussian model 

The distribution p( Gaussian ) ( J^) of J l3 for the Gaus- 
sian model is given in eq. (8). We set 



R - J ° 



(59) 



The thermodynamic quantity of the bond put be- 
tween the spin on node i and the spin on node j, 
[< b hond (i,j) >T ](^ ausslan ) j j Sj us j n g eqs _ (g, 10, 11, 
59) when /3 = /3 P , 

r , /. .n i (Gaussian) 

[< Obond(«,.?) >Tl R 

/oo poo 
■■■ n dj i™ n ^ (Gaussian) (^„) x 

■°° J -°°<lm> <lm> 



1 



E{Si} ^E<.->^- 



2 jv (2tt.7 2 ) j W 2 

/oo />oo 
■ ■ / fj rfJ; TO ^ PbondjSi, Sj,Jij) X 
-OO J —oo . r o l 



' <lm> {Si} 
J 2 

g- S<!m> YJT+PP T, < l m> JlmSlS 



= crf(/3 P J/V2) 



(60) 



where erf(x) is the error function of x. The thermody- 
namic quantity of the node degree for bonds at node i, 
[< fcbond(i) > T ](^ aussian ) ; i Sj us ing eqs. (8, 9, 10, 12, 59) 
when /? = Pp, 

[< fcbond(^) > T ](f a — ) = ert(p P J/V2) k(i) . (61) 
The thermodynamic quantity of the square of the node 



(Gaussian) 



is, 



degree for bonds at node i, [< k 2 ond (i) > t]r 
using eqs. (8, 9, 10, 13, 59) when /3 = P , 

\<r h 2 <i\ ^1 (Gaussian) _ 
[< fc bondW >T\ R - 

[eYi(f3pJ/V2)} 2 k(l)[k(l)-l} + 

erf(/3 P J/\/2)k(i). (62) 



We set 



^Gaussian) C£ p ) = erf (/3 P J/ ^2) . 



(63) 



Eqs. (60, 61, 62, 63) are formulated as eqs. (19, 20, 21), 
respectively. Therefore, the bond does not depend on the 
size of k(i). Using eqs. (14, 15, 23, 61, 62), we obtain 
the criterion of the percolation of the FK cluster for the 
Gaussian model as 



2crf(/3 P J/V2) > 2<fc> JV 



erf(/3 P J/y/2) + 1 <k 2 > N ' 



(64) 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 7 



Ineq. (64) is given by the inequality when the cluster is 
percolated. Ineq. (64) is given by the equality when the 
cluster is at the percolation transition point. 

From eqs. (18, 63), it is realized that there is the per- 
colation transition point for < f3p < oo. From eq. (64), 
it is realized that there is the percolation transition point 
for < a < 1. We approximate the error function erf(x) 
by 



erf(x) ~ y/l — cxp(— 4x 2 /tt) . 



Jo/ J is, using eqs. (59, 64, 65), 



Jo_ K (2-a) 2 
J V 2 4(1 - a) 



(65) 



(66) 



at the percolation transition point. The percolation tran- 



sition temperature T p 

rri(Gaussian) 



raussian 



is, using eqs. (59, 66), 
1 



k B . Itl i n (2-°)' 

2 111 4(l-a) 



(67) 



In Ref., 5 the internal energy on the Nishimori line for 
the Gaussian model, £;( Gaussian ) j i s derived. The internal 
energy ^(Gaussian) j s e q U i V alently 



^(Gaussian ) _ N Jo ^ 



>N ■ 



(68) 



The internal energy ^(Gaussian) us j n g e q S> (gg ? 53^ 



^(Gaussian) 




i, (2-a) 2 , 
2 ln 4(l-^ <fc>JV 



(69) 



We define the specific heat on the Nishimori line for 
the Gaussian model as (^(Gaussian) j n j^ e f 5 ^ ne U pp er 
bound of the specific heat c(Gaussian) j g (j er i vec j_ The up- 
per bound of the specific heat (7( Gaussian ) i s equivalently 



C (Ga 



< 



k B N(p P Jf 



< k> 



N 



(70) 



The upper bound of the specific heat c( Gaussian ) is, using 
eqs. (67, 70), 



^(Gaussian) < k B NlT ^ 



(2 



4(1 -a) 



< k > 



N 



(71) 



In Fig. 3, the percolation threshold of the FK cluster 
for the Gaussian model is shown. In (a) of Fig. 3, the 
relation between the aspect a of the network and Jo/ J is 
shown. Eq. (66) is used for showing (a) of Fig. 3. In (b) of 
Fig. 3, the relation between the aspect a of the network 
and the percolation transition temperature 7^ Gaussian ) j s 
shown. Eq. (67) is used for showing (b) of Fig. 3. In (b) 
of Fig. 3, J/k B = 1 is set. 

The complete graph is considerd as a case that a <~ 0. 
We set < k > N = N-l, < k 2 > N = (N-l) 2 , a = 2/(N- 
1), Jo — > Jo/N, and J — > J/y/W. From the settings, the 
model on the network becomes the SK model. 24 In what 
follows, our results for the SK model are shown. Jo/ J is, 
using eq. (66) for sufficiently large number of nodes, 



Jo 
J 



UN 



In 1 



iV 2 - AN + 3 



2N ™ 



10' 



10 l 



J /J 



10" 



10" 



10' 



10' 



10 l 



(a) 



(b) 



Un percolated 



Percolated 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 3. The percolation threshold of the FK cluster for the Gaus- 
sian model, (a) The relation between the aspect a of the network 
and Jo/ J is shown, (b) The relation between the aspect a of the 
network and the percolation transition temperature j'0 3ausslan ) 
is shown. J/fcs = 1 is set. 



at the percolation transition point. The percolation tran- 
sition temperature T p SK ^ is, using eq. (67) for sufficiently 
large number of nodes, 

t(SK) _ J 1 



'^ln(l+ ! 



N 2 -4,N+3 



J /27V 

The internal energy i?( SK ) is, using eq. (69) for suffi- 
ciently large number of nodes, 



(73) 



E (SK) 



J(iV-l) irN ( 1 

2 V~ ln V 1+ iV 2 -4iV + 3 



irN 



(74) 



For sufficiently large number of nodes, the right-hand 
side of eq. (71) reduces to 



k B N(N-l)n ( 1 

4 V N 2 -lN + 3 



k B TT 



■ (75) 



The upper bound of the specific heat C^ SK ^ is, using 



8 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



cqs. (71, 75), 



C (SK) < k BTT 

~ 4 



(76) 



The case for a = 1 is described. Using eq. (66), Jq/J = 
oo. Using eq. (67), Tp = 0. From eq. (16), it is realized 
that the network is at the percolation transition point. 
From Jq/ J = oo, the exchange interaction is only the 
ferromagnetic interaction. From Tp = 0, all the spins 
are parallel. Therefore, the cluster and the network are 
at the percolation transition point since the bonds are 
put on all the links of the network. 

In the result for the Gaussian model, an approximation 
formula of the error function, eq. (65), is used. In the 
result for the Gaussian model, it is necessary for the more 
precise estimation of the percolation threshold that the 
error function in eq. (64) is numerically estimated. 

7. Summary 

In the present article, the ±J model, the diluted ±J 
model, and the Gaussian model on the random graphs 
with artibary degree distributions are investigated. 

In the present article, for the bond that generates the 
FK cluster, [< b hond (i,j) >t]r, [< k hond (i) >t]r, [< 
fc bond(«) >t]r, [< k hond > t ]r, and [< fcg ond > t ]r are 
shown on the Nishimori line. They are rigorous even on 
a finite number of nodes. 

It is known that the internal energy, the upper bound 
of the specific heat, and so forth are exactly calculated 
on the Nishimori line without the dependence of the net- 
work (lattice). 5-9 In the present article, it is realized that, 
as a property for the Nishimori line, the bond does not 
depend on the size of the degree k{i) since the bias for the 
size of the degree k(i) does not appear in the statistical 
results of the bonds put between spins. 

In the present article, the percolation thresholds are 
shown. It is anticipated that ineq. (22) is true for suffi- 
ciently large number of nodes in the case that a condi- 
tion that the bond does not depend on the size of k(i) 
is satisfied. The condition is given in the present article. 
Ineq. (23) is derived from ineq. (22) and is the criterion of 
the percolation of cluster for spin models. Ineq. (23) gives 
the percolation threshold in the case that the condition is 
satisfied. The percolation thresholds are calculated based 
on ineq. (23). 

In the result for the Gaussian model, an approximation 
formula of the error function, eq. (65), is used. In the 
result for the Gaussian model, it is necessary for the more 
precise estimation of the percolation threshold that the 
error function in eq. (64) is numerically estimated. 

It is shown that the percolation transition temperature 
T p ±J ^ (eq. (32)) on the Nishimori line for the ± J Ising 
model on the present network coincides with the phase 
transition temperature j^ Ferro ) 26 > 27 for the ferromagnetic 



Ising model on the same network. For the ferromagnetic 
Ising model on the same network, the specific heat di- 
verges or jumps at j^ Ferro ) .26, 27 Q n ^ e o^er hand, from 
eq. (36), it is shown that the specfic heat neither diverge 
nor jump at Tp ±J ^ on the Nishimori line. 

The results for the infinit-range ± J model and the re- 
sults for the SK model are also shown. In Ref., 18 the 
percolation transition temperature of the FK cluster for 
the infinit-range ± J model is derived by using the ana- 
lytical solution of the SK model. The percolation transi- 
tion temperature of the FK cluster for the infinit-range 
± J model obtained in the present article agrees with the 
result for single replica case in Ref. 18 
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